Meshless acoustic analysis using a weakly singular Burton-Miller boundary integral formulation

The Burton-Miller boundary integral formulation is solved by a complex variable boundary element-free method (CVBEFM) for the boundary-only meshless analysis of acoustic problems with arbitrary wavenumbers. To regularize both strongly singular and hypersingular integrals and to avoid the computation of the solid angle and its normal derivative, a weakly singular Burton-Miller formulation is derived by considering the normal derivative of the solid angle and adopting the singularity subtraction procedures. To facilitate the implementation of the CVBEFM and the approximation of gradients of the boundary variables, a stabilized complex variable moving least-square approximation is selected in the meshless discretization procedure. The results show the accuracy and efficiency of the present CVBEFM and reveal that the method can produce satisfactory results for all wavenumbers, even for extremely large wavenumbers such as k = 10 000.     

高雷诺数下非混相Rayleigh-Taylor不稳定性的格子Boltzmann方法模拟

本文采用相场格子Boltzmann方法研究了竖直微通道内中等Atwoods数流体的单模Rayleigh-Taylor不稳定性问题,系统分析了雷诺数对相界面动力学行为以及扰动在各发展阶段演化规律的影响.数值结果表明高雷诺数条件下,不稳定性界面扰动的增长经历了四个不同的发展阶段,包括线性增长阶段、饱和速度阶段、重加速阶段及混沌混合阶段.在线性增长阶段,我们计算获得的气泡与尖钉振幅符合线性稳定性理论,并且线性增长率随着雷诺数的增加而增大.在第二个阶段,我们观察到气泡与尖钉将以恒定的速度增长,获得的尖钉饱和速度略高于Goncharov经典势能模型的解析解[Phys.Rev.Lett.200288134502],这归因于系统中产生了多个尺度的旋涡,而涡之间的相互作用促进了尖钉的增长.随着横向速度和纵向速度的差异扩大,气泡和尖钉界面演化诱导产生的Kelvin–Helmholtz不稳定性逐渐增强,从而流体混合区域出现许多不同层次的涡结构,加速了气泡与尖钉振幅的演化速度,并在演化后期阶段,导致界面发生多层次卷起、剧烈变形、混沌破裂等行为,最终形成了非常复杂的拓扑结构.此外,我们还统计了演化后期气泡与尖钉的无量纲加速度,发现气泡和尖钉的振幅在后期呈现二次增长规律,其增长率系数分别为0.045与0.233.而在低雷诺条件下,重流体在不稳定性后期以尖钉的形式向下运动而轻流体以气泡的形式向上升起.在整个演化过程中,界面变得足够光滑,气泡与尖钉在后期的演化速度接近于常数,未观察到后期的重加速与混沌混合阶段.     

非均匀波导中的声聚焦

理论研究了声波在非均匀波导中的空间聚焦问题,利用多模态导纳法构建波导内任意位置处声压与入射声压在模态域的映射关系,计算使声波聚焦于空间某位置时的最佳入射波,并画出了相应的聚焦声场.研究了三种非均匀情况:水平变截面波导、含散射体波导以及声速垂直变化波导.结果表明,当输入最佳入射波时,在非均匀波导中可以产生良好的单点或多点声聚焦效果,声波的聚焦过程充分地利用了波导结构及介质非均匀性对声波的散射作用.     

Positive ground state solutions for Schrödinger–Poisson system with critical nonlocal term in ℝ3

This paper is dedicated to studying the following Schrödinger–Poisson system $$\left\{\begin{array}{ll}-\mathrm{\Delta }u+V(x)u-K(x)\varphi |u{|}^{3}u=a(x)f(u),& x\in {ℝ}^3,\\ -\mathrm{\Delta }\varphi =K(x)|u{|}^5,& x\in {ℝ}^3.\end{array}\right.$$ Under some different assumptions on functions V(x), K(x), a(x) and f(u), by using the variational approach, we establish the existence of positive ground state solutions.     

80.5 MeV/u碳离子诱发铜靶的放射性剩余产物测量

高能重带电粒子能直接穿透靶原子核外电子层,与原子核发生直接碰撞,发生散裂反应,产生一系列具有放射性的剩余产物核.重带电粒子诱发靶材放射性剩余核与辐射防护和人员安全有着密切联系,当前,大部分剩余核产额主要依靠蒙特卡罗粒子输运程序进行模拟计算,其准确程度亟需通过实验测量进行准确评估.本文利用能量为80.5 MeV/u的(12)^C⁶⁺粒子对薄铜靶开展了辐照实验与伽玛射线测量,结合伽玛谱学分析方法,得出了辐照产生的18种放射性剩余产物的初始活度和产生截面值,并与PHITS模拟结果进行对比.结果表明,PHITS模拟程序对放射性剩余核种类的估计具有较高可靠性,在其绝对产额方面,与实验测量仍具有较大偏差.     

Flow physics of normal and abnormal bioprosthetic aortic valves

Flow physics of transvalvular flows in the aorta with bioprosthetic valves are investigated using computational modelling. For the efficient simulations of flow-structure-interaction in transvalvular flows, a simplified, reduced degree of freedom valve model is employed with a sharp interface immersed boundary based incompressible flow solver. Simulations are performed for normal as well as abnormal valves with reduced leaflet motion that models the effect of early leaflet thrombosis. The structure of the aortic jet and the hemodynamic stresses on the aortic wall are analysed to understand the hemodynamic impacts and possible long-term clinical implications of sub-clinical, reduced leaflet motion. The simulation results have shown that the reduced leaflet motion tilts the direction of aortic jet and generates stronger flow separation and re-attachment on the aortic wall downstream from the reduced motion leaflets. The modified flow pattern increases the wall pressure fluctuation and average wall shear stress on the downstream aortic wall, and results in the asymmetric oscillatory shear index distributions, which may have long-term clinical implications such as aortic wall damage and remodelling.     

Spectral estimates for saddle point matrices arising in weak constraint four‐dimensional variational data assimilation

We consider the large sparse symmetric linear systems of equations that arise in the solution of weak constraint four‐dimensional variational data assimilation, a method of high interest for numerical weather prediction. These systems can be written as saddle point systems with a 3 × 3 block structure but block eliminations can be performed to reduce them to saddle point systems with a 2 × 2 block structure, or further to symmetric positive definite systems. In this article, we analyse how sensitive the spectra of these matrices are to the number of observations of the underlying dynamical system. We also obtain bounds on the eigenvalues of the matrices. Numerical experiments are used to confirm the theoretical analysis and bounds.     

Arbitrage‐free call option surface construction using regression splines

In this work, we suggest a novel quadratic programming‐based algorithm to generate an arbitrage‐free call option surface. The empirical performance of the proposed method is evaluated using S&P 500 Index call options. Our results indicate that the proposed method provides a more precise fit to observed option prices than other alternative methodologies. Copyright © 2014 John Wiley & Sons, Ltd.     

砂岩方梁断裂过程区范围的实验确定方法

裂纹前端的断裂过程区是引起岩石非线性断裂及尺寸效应的主要原因。利用数字图像相关技术对砂岩开展了三点弯曲梁实验,获得观测区域高精度的全场位移和应变数据,根据断裂韧带区域水平位移和水平应变的分布特征,结合裂尖岩石颗粒变化的微观分析,提出采用裂纹尖端水平位移波动性和水平应变突变性所得到的波动系数和水平应变突变值,确定断裂过程区形状和临界尺寸的方法。结果表明:砂岩断裂过程区的形状为不规则的狭长带状区域,断裂过程区的临界长度为11~13mm,临界宽度为1.58~2.36mm。断裂过程区区域内形变在趋向裂尖时呈指数增加,但其单位区域内的形变增量呈波动状态。该方法能够更加准确判断岩石断裂过程区的范围,有助于分析岩石的非线性断裂特性。     

腔光子-自旋波量子耦合系统中各向异性奇异点的实验研究

通过耦合三维微波腔中光子和腔内钇铁石榴石单晶小球中的自旋波量子形成腔-自旋波量子的耦合系统,并通过精确调节系统参数在该实验系统中观测到各向异性奇异点.奇异点对应于非厄米系统中一种特殊状态,在奇异点处,耦合系统的本征值和本征矢均简并,并且往往伴随着非平庸的物理性质.以往大量研究主要集中在各向同性奇异点的范畴,它的特征是在系统参数空间中沿着不同参数坐标趋近该奇异点时具有相同的函数关系.在这篇文章中,主要介绍实验上在腔光子-自旋波量子耦合系统中通过调节系统的耦合强度和腔的耗散衰减系数两条趋近奇异点的路径而实现了各向异性奇异点,具体分别对应于在趋近奇异点时,本征值的虚部的变化与耦合强度和腔的衰减系数的变化会有线性和平方根不同的行为.各向异性奇异点的实现有助于基于腔光子-自旋波量子耦合系统的量子信息处理和精密探测器件的进一步研究.